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ABSTRACT 


The  spectral  chai*acteristics  of  random  modulated  waves  are 
investigated  with  particular  emphasis  on  the  determination  of  power 
density  spectra  of  sinusoids  that  are  frequency  modulated  by  band¬ 
pass  Gaussian  noise.  I’ox'Linent  parts  of  the  theory  of  random 
processes  necessary  for  spectral  analysis  are  reviewed,  and  a 
general  formulation  of  the  problem  of  determining  spectral  charac¬ 
teristics  of  random  modulated  waves  is  given.  The  principal  results 
of  this  investigation  are  given  in  Chapter  IV  where  the  variation  of 
FM  by  band-pass  noise  power  density  spectra  as  a  function  of  the 
modulation  parameters  is  treated  in  detail.  Limiting  cases  of  large 
and  small  modulation  indices  are  discussed,  and  a  technique  is 
developed  for  obtaining  useful  approximations  of  FM  by  band-pass 
noise  spectra.  These  results  are  verified  by  digital  compu’  ions 
and  experimental  measurements.  Finally,  spectra  resulting  from 
FM  by  the  combination  of  band-pass  noise  and  a  sinu.soid,  and  simul¬ 
taneous  AM  by  low-pass  noise  and  FM  by  band-pass  noise  are 
investigated. 


-.9- 


I.  INTRODUCTION 

Frequency  modulation  techniques  iiave  found  widespread 
application  in  electronic  communication  and,  more  recently,  radar 
and  radar  countermeasures  systems.  It  is  not  surprising,  there¬ 
fore,  that  considerable  effort  has  been  devoted  to  the  study  of 
frequency  modulation  fundamentals.  The  development  of  analytical 
techniques  for  determining  frequency  spectra  of  frequency  modulated 
waves  has  been  of  particular  interest,  where  the  modulating  wave¬ 
form  may  be  either  periodic  or  random. 

In  the  case  of  frequency  modulation  by  random  waveforms 
emphasis  has  been  placed  on  determining  spectra  resulting  from 
use  of  modulating  waves  that  are  obtained  by  filtering  white  Gaus¬ 
sian  noise  with  a  variety  of  low-pass  filters.  It  is  the  principal 
purpose  of  this  investigation  to  consider  modulations  which  consist 
of  Gaussian  noise  having  a  band-pass  characteristic;  in  particular 
the  case  pf  FM  by  white  Gaussian  noise  that  has  been  filter*"  '  by  a 
single  stage  L-R-C  filter  will  be  treated.  In  addition,  a  theory  for 
computing  spectra  resulting  from  FM  by  a  waveform  consisting  of 
a  sinusoid  and  band-pass  Gaussian  noise  combined,  and,  also,  FM 
by  band-pass  noise  with  simultaneous  AM  by  low-pass  noise  will 
be  given. 
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It  is  well  known,  and  almost  intuitively  obvious,  that  there 
exists  a  limiting  condition  wherein  the  spectral  shape  of  an  FM 
waveform  will  closely  resemble  the  probability  density  function 
of  the  modulating  voltage.  This  limiting  condition  is  associated 
with  an  increasing  modulation  index  (ratio  of  deviation  from  center 
frequency  to  deviation  rate).  Also  it  is  known  that  frequency  spectra 
of  amplitude  modulated  (AM)  and  FM  waves  are  almost  identical 
for  very  low  modulation  indices.  Therefore,  in  the  case  of  FM 
by  band-pass  noise,  an  interesting  transition  in  the  spectrum  occurs 
as  the  modulation  index  is  increased,  with  the  spectrum  going  Trom 
a  center  frequency  spike  with  two  sidebands  to  a  Gaussian  curve. 

It  will  be  possible  to  determine  the  modulation  index  reqviired  to 
produce  an  FM  wave  whose  spectrum  exhibits  an  essentially  Gaus¬ 
sian  shape  when  a  modulating  noise  filter  of  a  given  Q  is  employed. 

Formal  analytical  methods  of  computing  FM  spectra,  in  all 
but  the  cases  of  relatively  simple  modulating  waveforms,  invari¬ 
ably  lead  to  complicated  macbematical  expressions  which  are 
exceedingly  difficult  to  evaluate  explicitly,  and  give  little  insight 
into  the  interrelationship  of  the  various  modulation  parameters., 

Thus  one  is  forced  to  seek  approximate  expressions  which  are 
capable  of  giving  reasonably  accurate  results  for  the  parameter 
range  of  interest.  This  state  of  affairs  is  not  unexpect'jd  when 
dealing  with  frequency  modulation,  but  is  clearly  a  result  of  ti>e 
fact  that  frequency  modulation  is  a  nonlinear  process.  In 
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noiilinear  analysis  approximate  solutions  are  generally  the  only 
useful  solutions  that  can  be  obtained.  Of  course,  in  order  to  have 
confidence  that  approximate  solutions  will  provide  useful  results, 
it  is  necessary  to  determine  the  conditions  (or,  equivalently,  the 
range  of  parameters)  under  which  the  approximations  are  valid; 
and  also,  insofar  as  possible,  approximations  should  be  verified 
by  comparison  with  experimental  results  obtained  from  the  actual 
or  simulated  nonlinear  system  and  by  evaluation  of  exact  calcula¬ 
tions  where  the  parameters  have  been  fixed  to  make  computation 
possible.  In  this  investigation  of  FM  spectra,  approximations  will 
be  suppor^d  by  both  experimentally  determined  and  digitally 
computed  spectral  data; 
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II.  A  SUMMARY  OF  THE  THEORY  OF  RANDOM  PROCESSES 

The  detei'mination  of  epcctral  characteristics  of  random 
modulated  waves  is  essentially  an  applied  problem  in  the  theory 
of  random  processes.  Accordingly,  a  summary  of  applicable 
portions. of  the  theory  will  be  given  in  order  to  provide  a  basis  for 
a  detailed  formulation  of  the  power  spectrum  problem. 

Extensive  literature  exists  on  the  general  subject  of  random 
(or  stochastic)  processes  with  books  by  Doob  (1)^,  and  Davenport 
and  Root  (2)  being  notable  examples.  The  former  is  a  rigorous, 
inathematical  text,  while  the  latter  contains  an  applied,  engineering 
treatment  of  the  subject.  Throughout  the  literature  there  is  con¬ 
siderable  variation  in  the  terminology  and  notation  that  are  employed 
to  denote  quantities  and  relationships  that  are  essential  for  power 
spectrum  computations.  The  notation  introduced  in  this  chapter 
will  be  employed  throughout  the  remainder  of  this. report,  and 
defined  quantities  will  be  underlined  for  emphasis.- 

A.  THE  RANDOM  PROCESS 

A  random  process  can  be  simply  defined  as  a  collection 
(or  ensemble)  of  random  time  functions.  A  member  function  of  the 
ensemble  can  be  written  as  x^(t)  where  t  is  a  continuous  parameter. 


Numbers  in  parentheses  refer  to  references  given  on  page  84. 
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denoting  time,  which  ranges  from  -oo  to  +oo,  and  a  is  a  typical 
point  in  a  probabiUty  measure  space  Si.  A  classic  example  of  a 
random  time  function  that  is  of  central  importance  in  electronics 
is  the  noise  voltage  produced  by  a  resistor  due  to  thermal  agitation 
of  electrons.  A  collection  of  the  noise  voltages  produced  by  all 
possible  physical  realizations  of  a  resistor  of  R  ohms,  for  example, 
would  comprise  a  random  process.  Since  it  is  not  feasible  to 
specify  each  time  functio..  of  a  random  process  exactly,  various 
averages  rriust  be  defined  in  order  to  describe  the  properties  that 
can  be  expected  to  apply  to  an  arbitrary  member  fiinction  of  the 
ensemble.  This  will  lead  the  way  toward  computing  some  important 
properties  of  the  output  of  a  sy.stem  which  has  a  random  input.  Tb  > 
"system"  that  will  be  considered  in  this  report  is  an  oscillator 
which  has  the  capability  of  being  modulated  in  both  amplitude  and 
phase,  and  we  shall  be  interested  in  studying  spectral  properties 
of  the  output  voltage  for  inputs  which  consist  of  various  random 
modulations . 

1 .  Probability  Distribution  and  Density  Functions 

The  random  time  functions  x  (t)  evaluated  at  a 

a 

specific  time,  t.  ,  comprise  a  random  variable  which  is  defined 
over  the  space  n.  (Note:  Usage  dictates,  the  use  of  the  word 
variable  instead  of  function  as  would  be  more  appropriate.)  The 
n^^’  order  cumul.ativc  probability  distribution  function  of  tlu-  set 
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o£  random  variables  x  (t, ),  x  (t-),  •  •  • ,  x  (t  )  is  d-moted  P  ,  and 

Cl  X  Cl  ^  01  n  iix 

defined  as 


^’nx^’^l*  ’‘Z*  ■  ■  ■  *  ^n‘  *1'  *2’  ‘  "  ’  ‘^n^ 

=  Prob  [{a|x^(t^Ux^.  ....  x^(gs  x^j  J 

where  the  second  expression  should  be  read  as  the  probability  that 

x^(t)  evaluated  at  time  t^  will  be  less  than  or  equal  to  the  value 

x^  and  simultaneously  x^(t2)  is  less  than  or  equal  to  x^t  and  so 

on  to  gtjj)  s  x^.  An  equivalent  interpretation  is  that  the  set  of 

points  o  which  correspond  to  those  member  functions  of  the 

ensemble  that  meet  the  condition  j^x^(tj)  s  x^,  •  •  •  ^  J 

has  a  measure  given  by  P  .  The  function  P  ,  has  the  usual 

nx  ni 

properties  associated  with  probability:  for  example,  0  s  P^^^  s  1 
and  P^^  is  monotonic  nondecreasing.  In  spectral  analysis  the 
second  order  probability  distribution  function,  P2x»  will -be  of 
central  interest. 

The  probability  density  function,  is  defined 

as  the  mixed  partial  derivative  of  P^^^^  when  this  derivative  exists. 
Thus 


•t  )=, 
n'  * 


Xnilf.g 
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The  joint  cumulative  probability  distribution  function  of  order  (nx, 


juy)  for  the  random  processes  x^(t)  and  .\s  given  by 


=  Prob  [  {a|  x^(t^)^  x^.  .  •  •x^(y  ^  x^;  y^(r^)^  y^.  ■  '  -yjr^jsy^}] 


Two  random  processes  are  called  independent  to  order  (nx,  my)  if 


Pnx,my<’'l*^'-yn'  yi*"‘*ym'  ‘l*"‘*‘n=  ’’l'"-  V 


=  P^y(yi.---y^!  Tj.---  T^) 


or,  in  words,  the  joint  distribution  of  x  and  y  is  equal  to  the 
product  of  the  individual  distributions. 

In  view  of  the  fact  that  the  averages  that  are  important 
in  spectral  theory  can  generally  be  obtained  itsing  first  and  s  -ond 
order  distributions,  henceforth  higher  ordered  distributions  will 
not  be  employed  in  the  development  to  follow.  This  does  not 
represent  any  real  loss  in  generality,  however,  since  the  extension 
to  higher  order  statistics  is  straightforward. 

One  last  definition  pertaining  to  distribution  functions 
and  al.'0  density  functions  is  the  notion  of  stafcionarity ,  A  distribu¬ 
tion  is  said  to  be  stationary  to  second  order  if  the  distribution 
function,  P2^»  corresponding  to  x^(t),  is  equal  to  the  distrih  .ion 


corresponding  to  x^(t+T)  where  t  is  a  fixed  parameter.  That  is, 


2x'  l'"2*  -1 


This  reduces  the  number  of  variables  by  one  since  r  can  be  taken 

et.  .al  to  -t^  and  thus  is  a  function  of  x^  ,  x^,  and  1^ .  A 

tVi 

similar  definition  applies  for  n  order  distributions.  Also  first 
order  stationary  distrJbuti..'.ia  are  independent  of  time . 


B.  STATISTICAL  AVERAGES 

Since  the  random  process  is  characterized  as  a  function 
of  two  variables,  it  is  apparent  that  two  distinct  types  of  averages 
are  possible.  These  are:  (1)  averages  over  the  ensemble  variable, 
a,  at  a  fixed  time,  t^;  and  (2)  averages  over  time  for  a  given 
member  function  of  the  ensemble. 


1,  Ensemble  Averages 


The  ensemble  average,  or  expected  value,  of  the 


random  variable  x  is  defined  as 


=  I  V  Plx 

r» 

where  E  is  the  expected  value  operator.  Various  moments  of  the 
random  variable  x  arc  given  by  E(x”),  Of  principle  importance 
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Finally  it  should  be  noted  that  the  expected  value 
inie=i-al  reduces  from  a  Stieltjes  integral  to  an  ordinary  Riemann 
integral  if  has  a  continuous  derivative.  That  isi 


E[f(x^)]=  J  J  f{Xi)Wi^(x^,t^)dx 

J  •'-CD  “-00 


Jump  discontinuities  in  the  monotonically  nondecreasing  distribution 
function  can  usually  be  accounted  for  by  using  delta  functions 

with  the  density  function,  If  g  is  a  function  of  the  two  vari¬ 

ables  x^  and  then  the  average  of  g  becomes 


.00  .,00 


E|^g(x^,X2)J  =  J  J  gix^.x^)  W2j^(Xj,  X2J  t^,t2)  dx^  dx2 


-CO  -00 


Note  that,  in  general  E|f(x^)J  is  a  function  of  t^  and  E^g(x^,X2)J 
is  a  function  of  t,  and  t,;  however,  if  the  distributions  W.  and 

14  ,1 X 

Wj,^  are  stationary,  E|f{x^)  J  is  time  independent  and  E|g(x^,x^)J 
depends  only  on  the  time  difference  t^  -  t, , 


An  average  of  considerable  importance  in  probability 
theory  is  the  characteristic  function,  which,  for  the  second  order 
case,  is  defined  as 


E 


[_exp 


(jv 


I’^l 


+  jv 


z’^z 


>]  = 


t,,y 
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If  W«  existf?  then  the  characteristic  function  becomes 
2x 


CO  00 

Mzx<''i-''2=4'V=  J  J  expOv^x^+jv^Xj,)  W^^{Xj,X2:t^,t2)  dx^ 

-00  -00 

Thus  the  characteristic  function  is  the  Fourier  transform  of  the 
density  function  {except  for  the  sign  of  the  exponent  of  e). 

The  characteristic  function  has  great  utility  in 
computing  the  density  function  of  the  sum  of  two  random  variables 
when  the  individual  density  functions  are- known.  Consider  the 
first  order  case  of  the  sum  of  x^  and  y^  v/ith  and 

W^y(yi)  given.  The  subscripts  wili  be  dropped  for  convenience 
(a  standard  practice)  and  stalionarity  of  and  will  be 

assumed.  Let  z  =  x  +  y  and  the  problem  is  to  find  The 

characteristic  function  of  z  is 

^Iz  ~  (jvz)J  =  Ej^exp  (jvx  +  jvy)J  . 


If  X  and  y  are  indcpendenl  tiien 


M,  =  M,  M, 
Iz  lx  ly 


The  density  function  of  z,  W.  can  now  be  computed  from 

y  2 


A 

~  Tv  I  M.  (v)  exp  (-jvz)  dv  . 
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2.  Time  Averages 

The  average  over  time  that  is  performed  on  members 
of  the  ensemble  is  defined  as 


where  the  time  average  operator  is  denoted  by  A.  The  two  integrads 
usually  reduce  to  the  single  '-.legral  formula 


'  Ti»  rr  J  \'‘i  ■“ 

-T 

in  view  of  the  fact  that,  when  the  two  integral  expression  exists, 
the  single  integral  expression  also  exists  and  they  are  equal.  It 
is  possible  that  the  single  integral  expression  may  exist  when  the 
two  integral  expression  does  not  for  certain  functions,  and  in  order 
to  obtain  the  desired  property  of  time  invariance  for  the  linear 
operator  A,  the  two  integral  formulation  of  time  average  must  be 
used. 

The  time  average  of  a  member  function  x^(t)  of  a 
random  process  will,  in  general,  depend  on  the  ensemble  variable 
a.  When  the  time  averages  of  the  member  functions  are  independent 
of  a,  the  random  process  is  said  to  be  uniform.  Actually  the 
character  of  the  integral  used  for  ensemble  averaging  is  such  that 
the  results  obtained  are  not  changed  if  at  most  a  ,  ountable  infin’ty 
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of  member  functions  do  not  have  the  property  of  uniformity.  If 
the  notion  of  a  random  process  is  to  be  used  in  tiic  formulation  of 
a  mathematical  model  for  a  physical  system,  the  assumption  of 
uniformity  must  generally  be  made,  in  that  only  the  member  function 
of  the  random  process  is  e.vailable  for  experimental  observation. 

3.  Equivalence  of  Time  and  Ensemble  Averages 

Suppose  ensemble  and  time  averaging  are  performed 
on  the  function  gCx^.x^)  with  reversed  order  as  follows: 

AE  j^g{x^,  x^lj  =  EAj’g(x^,X2j  . 

Various  necessary  conditions  for  the  existence  of  the  integral 
A^g(x^,  X2J  are  given  by  theorems  on  ergodic  theory.  Now  if  the 
process  is  uniform,  the  relation  reduces  to 


AE[g(x^,x,)‘j  =  Aj^g(x^,X2)J 

since  A[g(x  ie  independe:it  of  the  ensemble  variable  and, 

therefore,  the  ensemble  average  leaves  A|^g(Xj,X2)J  unchanged. 
Also  if  stationary,  then 


E  [g(Xj.  x^)]  =  A|*g(x^,  X2,)j 


-22- 


This  relationship  provides  th^|jnk  between  ensemble  averages 
computed  from  the  random  process  model  of  a  physical  system, 
and  the  measurements  that  may  be  made  on  an  actual  system 
using  time  averages.  An  excellent  and  thorough  treatment  of  the 
general  subject  of  the  relationships  between  ensemble  and  time 
averages,  which  is  reviewed  only  briefly  here,  may  be  found  in 
a  paper  by  W.  M.  Brown  (3). 

C.  CORRELATION  FUNCTIONS  AND  POV/ER  SPECTRA 

With  the  material  covered  in  prececding  sections  of  this 
chapter  available,  a  theory  for  spectral  characteristics  of  random 
processes  can  be  given  without  difficulty.  The  theory  presented 
here  will  be  based  on  ensemble  averages;  however,  equivalence 
with  the  spectral  theory  obtained  for  time  averages  may  be  qiUckly 
demonstrated  with  the  aid  of  the  results  of  the  previous  section. 

1 .  Ensemble  Cov**elation  Functions 

The  autocorrelation  function  for  x  (t)  over  the 
- — - -  a 

ensemble  is  defined  as 

Rx(M«4)  =  e[x^’’‘z]  =  J  J  V  ‘^’‘2  • 

-00  -00 
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If  is  stationary,  then  is  a  function  of  the  time  difference 

t,-tj  which  will  be  denoted  as  t-  If  a  random  process  has  the 

property  that  its  correlation  function  depends  only  on  the  time 

t2T  — 7 

=  x^  is  independent  of  time;  then 
the  process  is  called  stationary  in  the  wide  sense.  Cross¬ 
correlation  functions  are  defined  for  two  random  processes  x^(t) 
and  y^(t)  as 

however,  we  shall  not  have  need  of  functions  of  this  type;  conse¬ 
quently  the  terms  autocorrelation  and  correlation  will  be  synonymous 
for  the  remainder  of  this  report. 

An  important  property  of  correlation  functions  can 
bo  easily  derived  by  observing  that  the  expected  value  of  a  non- 
ncgativc  quantity  is  non-negative .  Thus 
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thus 

and  for  stationary  random  processes 

R^{0)  s 

since 


2 .  Power  Density  Spectra  of  Random  Processes 

The  power  density  spectrum,  Sj^(w)>  of  a  random 
process  x^(t)  is  defined  as  the  Fourier  transform  of  the  correla¬ 
tion  function  R  (t)  where  the  transform  exists.  Note  that  x  (t) 

X  a 

is  assumed  to  be  wide-sense  stationary.  Since  a  function  is  re¬ 
quired  to  be  absolutely  irtegrable  to  insure  existence  of  its 
Fourier  transform,  unless  R{t)  goes  to  zero  as  t  goes  to  infinity, 
the  power  density  spectrum  will  not  exist  in  the  strict  sense.  Then 
the  correlation  function  is  related  to  a  cumulative  power  spectrum 
S^(w)  by 

4  pOO 

~  Tn  1  ■ 

-CO 
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This  last  result,  which  is  essentially  the  famous  Wiener-Khinchin 
Theorem,  provides  the  basic  procedure  for  coniputing  power 
spectra  of  random  waveforms.  As  in  the  case  of  probability  dis¬ 
tribution  and  density  functions,  it  is  usually  possible  to  introduce 
delta  functions  to  account  for  jump  discontinuities  in  the  cumulative 
power  spectrum.  The  .appearance  of  delta  functions  in  power  density 
spectra  accounts  for  power  corresponding  to  discrete  frequencies 
including  w  =  0  {or  dc).  it  should  be  noted  that  there  exist  correla¬ 
tion  functions  for  which  the  cumulative  power  spectrum  has  the 
property  that  it  is  continuous  everywiiere  except  at  a  finite  set  of 
points;  however,  its  derivative  there  is  zero  (the  so  called  singular 
case).  In  this  situation,  the  delta  function  cannot  be  used;  however, 
such  cases  are,  fortunately,  not  commonly  encountered  in  engineer¬ 
ing  practice.  The  use  of  power  density  spectra  containing  delta 
functions  will  be  adequate  for  the  power  density  spectra  that  will 
be  computed  in  this  report. 

It  remainr  to  give  a  Justification  for  the  name  power 

density  spectrum  for  the  fu.uction  which  has  been  rather  arbitrarily 

defined  as  the  Fourier  transform  of  an  ensemble  correlation  function. 

If  the  inverse  Fourier  transform  is  used  to  express  in  terms 

of  S  (w)  then 
x'  ' 

,  ,oo- 

"  Zir  J  ®x^“^ 

-00 
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The  correlation  function  evaluated  at  t  =  0  is  simply  the  mean 
square  value  of  fhus 

-00 

Also  it  can  be  shown  (2)  that  the  correlation  function  of  the  output, 
y^(t)  of  a  linear  time  invariant  system  can  be  given  in  terms  of 
the  correlation  function  of input,  by  the  triple  convolution 

Ry(T)  =  h(t)  *  h^(t)  *  R^(t)  , 

where  h(t)  is  the  system  impulse  response  and  h^(t)  is  the  complex 
conjugate  of  h(-t).  In  the  frequency  domain  this  becomes 

Sy{<*>)  =  S^(U>)  , 

where  H(j.w)  is  the  system  transfer  function  (Fourier  trcuisform  of 
h(t)). 

If  the  system  characterized  by  .H(jw)  is  assumed  to 
be  a  narrow  band  filter  with  unit  response  for  s  u 
aero  elsewhere  then 

“2 

“y"”  J  d„  . 
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Now  since  y  is  always  positive  and  since  and  W2  can  be 
closely  spaced  at  any  point  along  the  frequency  axis  w,  it  follows 
that  S^(w)  mast  be  a  non-negative  function  of  w.  Also  it  can  be 
seen  that  S^u>)  is  an  even  function  by  observing  that 


(*°°  1*°° 

SJu)  =  J  R^{r)  exp  (-Jwt)  dr  =  J  R^{r)  cob  ut  dr 


“00 


-00 


since  R^ir)  is  a  real,  even  function;  and,  therefore. 


,.00  00 

=  J  R^if)  cos  (-a>T)  dr  =  J  cos  wt  dr 


-00 


-00 


Thus,  S^(w)  is  seen  to  correspond  to  the  distribution  of  power  with 
frequency; 

A  sinhilar  theory  can  be  developed  for  time  averages 
by  defining  the  autocorrelation  function  over  time  as(5^(T)  = 

A  P x(t)  x(t+r)J  for  uniform  random  processes.  The  time  power 
density  spectrum  is 


"  J  d  ir)  exp(-JcoT) 


dr 


If  x^(t)  is  not  stationary  in  the  wide  sense  then  R^(t,T)  = 
E|x^(t)  x^(t+T)J  can  be  averaged  over  time  to  eliminate  the  t 
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dependence  and  a  power  spectrum  is  computed  as 


.00 


=  J  exp  (-jwT)  AE  l^x^(t)  x^(t+T)J  dr 

-CO 

This  power  density  spectrum  is  also  denoted  since 


AE^x^tt)  x^(t+T)J  =  A  jx^(t)  x^(t+T)J 


Now  if  x^(t)  is  wide  sense  stationary  then 


E  [x^(t)  xJt+T)J  =  A  [x^(t)  X^(t+T)J 


and 


=  S  (w) 


X'  '  X 


D.  THE  GAUSSIAN  RANDOM  PROCESS 

The  Gaussian  random  process  is  of  particular  importan 
in  the  analysis  of  physical  systems  because  of  its  frequent 
occurrence.  For  example,  many  noise  voltages  encountered  in 
electronic  systems  are  characterized  by  either  Gaussian  distribu¬ 
tions  or  nonlinear  transformations  of  Gaussian  distributions.  The 
wide-spread  occurrence  of  phenomena  that  have  Gaussian  properties 
is  predicted  by  the  central  limit  theorem  which  states  ths.'c  distribu¬ 
tion  of'the  sum  of  n  independent  random  variables  becomes  Gaussian 
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as  n  tends  to  infinity.  This  would  serve  to  explain  the  fact  that 
the  noise  voltage  of  an  ohmic  resistor  is  Gaussian  in  that  it  is 
made  up  of  the  combination  of  a  large  number  of  impulses  due  to 
thermally  ihduced  motion  of  individual  electrons. 

The  first  and  second  order  Gaussian  density  functions 
and  the  corresponding  characteristic  functions  will  be  recorded 
here  for  later  use.  It  will  be  assumed  that  the  random  process 

has  a  zero  mean  valuCt  is  E|^x^(t^)J  =  0 


exp(-— ^r)  =  exp  (- -i-TT-i- ) 


2  2 
<r 

T 


exp 


2  2  ,  ,  2  2 
Xj  -  ZR^{t^,t^)  x^x^  +  Xj> 


]}■ 
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An  important  property  of  Gaussian  distributed  quantities 
is  the  fact  that  if  the  input  to  a  linear  system  is  Gaussian,  then  the 
output  is  Gaussian.  In  other  words  linear  transformations  of 
Gaussian  noise  give  Gaussian  noise.  (See  Reference  4  for  details. ) 
It  is  true,  however,  that  a  linear  transformation, may  well  result  in 
Gaussian  distributions  that  are  non- stationary. 
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m.  FORMULATION  OF  THE  RANDOM  MODULATED  WAVE 
PROBLEM 

A.  THE  ENSEMBLE  CORRELATION  FUNCTION  OF  A 

GENERAL  RANDOM  MODULATED  WAVE 

A  general  expression  for  the  random  time  function  whose 
spectral  properties  wHl  be  investigated  is  as  follows; 

v(t)  =  C(t)  cos  j^o)^t  +  4)(t)  +  0^  , 

Here,  specifically,  v(t)  represents  the  instantaneous  voltage  out¬ 
put  of  a  sinusoidal  oscillator,  of  radian  frequency  whose 
amplitude  is  modulated  by  C(t)  and  whose  instantaneous  radian 
frequency  is  given  by  +  ^{t)  with  0  being  a  fixed  arbitary 
phase  angle.  Now  if  C(t)  and  <})(t)  are  random  variables,  then 
v{t)  is  a  random  variable.  The  ensemble  parameter,  a,  will  not 
be  indicated  explicitly;  however,  it  is  to  understood  that  v(t), 

C{t),  and  <|>{t)  represent  member  functions  of  random  processes. 

The  correlation  function  of -v(t)  will  now  be  formulated; 
however,  the  question  ae  lo  '.''idc  sense  stationarity  of  v(t)  will  be 
left  open  for  the  moment,  and  R^  will  be  indicated  as  a  function 
of  both  t^  and  or  equivalently  t  =  t^  and  T=  t^-t^. 
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Ry(t,  T)  =  E  ^V{t)  V(t+T)j 

=  E|c(t)  C(t+T)  COP  j^w^t+4i{t)+6j  cos  |^a>^t+u^T+<t>(t+T)  +  0  J I 

Using  the  trigonometric  identity  for  the  product  of  two  cosine 
functions  along  with  the  linearity  property  of  the  expected  value 
operator  E  gives 

R^(tiT)  =  e|^  C(t)  C(t+T)  cos  [Zw^t  +  w^T  +6(t+T)  +  (i>{t)  +  20  J  | 

+  e|-j  C(t)  C(t+T)  cos  +  <|>(t+T)  -  iMt)]  I  . 


The  arbitrary  phase  angle  0  is  assumed  to  be  uniformly  d^tribu^d 
over  the  interval  0  s  0  s  2ir,  that  is 


1 


Os  0  s  2ir 
for  all  other  0 


The  required  joint  probability  density  function  for  the  expec^d 
value  integral  would  have  the  form  ^2+  ^10  assump¬ 

tion  that  the  random  variables  Ct  if*  ^nd  0  are  independent.  In 
carrying  out  the  indicated  expected  value  operations  to  obtain  the 
correlation  function  R^{t,T)*  the  term  which  involves  the  phase 
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angle  0  will  vaniah.  This  can  be  seen  by  noting  that  the  expression 


e|.j  C(t)  C(t+T)  cos  j^Zw^t  +  +  <t'(t+r)  +  <Mt)  +  2eJ  j. 


takes  the  form 


^3(1,7)  J^’'cos[p(t.T)4  20]  W^0(G)d0 


-00  -  00  -00  -00 


and 


dC|  di|)^  d4)2 


A  »2it 

r  cos  (P  +  20)  dO  =  0 
0 


Thus  the  correlation  function  of  v(t)  reduces  to 

R^(t,T)  =  C(t)  C(t+T)  cos  [wj  +  <|.(l+T)  -  <|.{t)]  j  . 

In  order  to  proceed  further  the  functions  C(t)  and  i))(t)  must  e 
specified.  The  following  cases  will  be  considered  in  detail  in  the 
chapters  to  follow; 

(1)  Case  I.  FM  by  band-pass  noise  only.  C(t)  is  a 
constant  C^,  and  <Mt)  is  obtained  by  integrating  band-pass  noise, 


Mt)  =  D  V  (f)  df  , 

-CO 


that  is 
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I 

> 

I 

» 

I 

I  °  «en.lUvity  of  the  einueoldal  o.olUalor 

I  in  Mdian./volt  and  v„(t)  i.  a  atationany  Gaua.ian  noiee  voltage 

,  With  zero  mean. 

(2)  Case  H,  FM  by  band-pass  noise  combined  with  a 
sine  wave.  This  is  the  same  as  Case  I  except  that  ^t)  is  now 

^  ^  J  (V'  +  . 

-CO 

where  and  o,^  are  the  peak  value  and  radian  frequency  of  a 
sme  wave  respectively,  and  o  is  a  uniformly  distributed  random 
phase  angle.  The  band-pass  noise  spectrum  will  be  centered  at 


(3)  Case  III.  FM  by  band-pass  noise  with  simultaneous 
am  by  Io«-p...  noiho.  The  function  +|l)  i„  ,h„  ..  c,.e  f 

but  C|t)  i,  no*  C(t)  =  c„  +  v^(t)  where  v,|t)  i.  a  et.tionary 
Caueeten  noiee  voltage  with  aero  mean  and  having  a  low-pae,  power 
deneity  epectrnm.  Ai.o  v^(t)  end  v^(t)  a.eumed  to  be 
indenendent. 
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IV.  FM  BY  BAND-PASS  NOISE 

The  development  of  this  chapter  for  the  case  of  FM  by 
band-pass  noise  represents  the  principal  contribution  of  this 
report.  A  method  of  approximation  with  which  useful  informa¬ 
tion  concerning  the  shape  of  the  power  density  spectrum  as  a 
function  of  the  essential  modulation  parameters  will  be  explored 
in  detail. 


A.  SPECIFICATION  OF  THE  CORRELATION  FUNCTION 
IN  TERMS  OF  THE  MODULATING  NOISE  SPECTRUM 

The  correlation  function  for  FM  by  band-pass  noise 
was  derived  in  Chapter  Ill  as 

C  ^ 

Rv(t,T)  =  -5"  E  {cos  [u,^t+  <Ht+T)  -  +{1)]  }  , 

with 

iMt)  =  D  J  v^(t')  df  . 

-00 

If  the  power  density  spectri-m  S  {<*>),  of  the  modulating  noise  is 

n 

specified,  then- the.  cor  relation  function  can  be  evaluated  as  will 
be  shown  subsequently.  Also  the  issue  as  to  the  wide  sense 
stationarity  of  v(t),  or  equivalently  the  question  of  time  indepen- 
dependenco  of  R^,  will  be  settled. 
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The  correlation  functioni  Ry(ti  t),  can  be  written  in  the 
following  form  where  the  equivalent  time  variables,  t^  =  t  and 
-  t+T,  have  been  used  in  place  of  t  and  t,  so  that  (Mt+r)  = 

^2  **“1  "Wti)  =  <Mt)- 

C  ^ 

Ry{t,T)  =  Re  |exp(Ju^T)  E  j^exp  (j<|>2  "  j'I'i)]]  • 

The  symbol  'Re  denotes  "real  part  of".  From  the  definition  of  the 
second  order  characteristic  Junction  of  this  is  seen  to  be 

C  2 

Since  <}>  is  a  linear  transformation  of  a  Gaussian  process,  it  too 
will  be  Gaussian,  but  not  necessarily  stationary;  therefore 

— z  — Z  ^ 

The  terms  4)^  ,  i})^  ,  and  R^(t|»  t^)  =  9^"?^  will  now  be  computei 
from  the  defining  relation  for  i^,  that  is 

pt 

«Mt)  =  D  J  v^(t')dt'  . 

-oo 


Lat  y-x  =  u  and  y  =  v.  This  change  of  variables  maps  the  x,  y 
plane  into  the  u,  v  plane  as  shown  in  Figure  1. 
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FIGURE  1  KIAPPING  OF  x,  y  TO  u,v  PLANE 


The  integral  for  R^(t,T)  is  obtained  for  T  ^  oo.  In  terms  of  the 
new  variables  ii  and  v,  the  integral  becomes 


u+t 


=  D 


u+t. 


'  0  r  1  1  T  r  1  I 

■  J  (U)  1^  j  dv  du  +  J  (u)  J  dv  J 


du 


-39- 


Similarly  4)^  and  become 


r  (*0 

U+tj 

{J  V'“’ 

J  dv 

and 


+7  = 


-00  -  -00 


u+t. 


-00 


'{f  ^v]au}. 

-00  -CO  0  -oo 

<j>^  +  (|>2  -  2  4)^<j)2  J  can  now  be  evaluated. 


♦j  +4>->  ffl  r  1  r  1 

J-^=fR  (m[j  dvldu4j  R  (u)rj  dvjdu 

^  rOO  ^  -00  "^0  *■  -00 

0  r  -1  00  r  ^2  T 

+  J  Ry  (u)  J  dv  du  +  J  (u)  j  dv  du 

-oo  ”  -oo  0  -oo 


r 

-u+t. 

1%  * 

-U«J 

1 

=  J 

J  dv  + 

J  “■'  + 

J 

0  "  ^ 

-00 

-00 

-00 

J 

-00 

where  use  has  been  made  of  the  fact  that  R  (u)  is  an  even  function. 

n 
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At  this  point  the  correlation  function  R  (u)  is  expressed  in  terms 

'^n 

of  the  corresponding  power  density  spectrum  (w)  to  get 

n 


+  00  T 

"-j - ^ —  =  —  J  (w)  j  (t-u)  cos  wu  du  dw  , 

-00  ”  0 


D- 


where  the  Fourier  cosine  transform  expression  for  S  (w)  is  used, 

n 

that  is 


1 

"Iff  wr  dw 

n  "  n 


-oo 


The  integral  over  u  is  easily  evaluated  to  get  finally 


dw 


•00 


Thus  the  correlation  function  for  v{t)  becomes 
^  2 


.os  w^T  exp  [-  r  (w) 

-00  ” 


R^(t)  =  — j—  cos 


Note  that  R^  has-been  shown  to  be  independent  of  time. 
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B.  THE  POWER  DENSITY  ISPECTRUM  OF  WHITE  NOISE 
PROCESSED  BY  A  SINGLE  STAGE  LRC  FILTER 

The  modulating  voltage  v^(t)  is  obtained  by  filtering 
"white"  Gaussian  noise  by  a  single  stage  LRC  filter  such  as  that 
shown  in  Figure  2  below. 


FIGURE  2  SINGLE  STAGE  LRC  FILTER 


The  transfer  function  of  this  filter  is 


H(8) 


’m 


"C" 


8 

T“ 

Rir 


where  g„  is  the  transconductance  of  the  tube.  Let 
m 


1 

“i  =  rC 


and 
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The  amplitude  frequency  response  squared  becomes 


!H{ju.)| 


u  2  2 

Hq  0, 


,  2  2.“^  2  2 
(Wq  -  w  )  +  «^  u 


with  H, 


2  8 


m 


If  the  power  density  spectrum  of  v.  (t)  is  white,  that  is,  S  (w)  =  N 

in 

for  all  «,  then  the  power  density  spectrum  of  the  output  v^(t)  1* 


Sv  <“)  =  — ^ - z — ■ — 

''n  ,  2  2.^^  2  2 

(“O  -  “  )  + 


The  mean  square  value  of  v^(t)  is  given  by 


N  .00 


■ — -y-  j  w  It  w  u 

%  =  4  I  ^  =  -?^r-  J  ;  2-  z/,  2--  g 

n  *’  (Un  -  W  )  +  W,  U) 

-00  -oo  '0  '1 


du 


Evaluation  of  the  integral^  yields  ^  . 


Thus 


-TT,  _  "o'" 


3 


See  Appendix  A. 


FIGURE.  BAND-PASS  MODULATING  NOISE  SPECTRA 
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Let  m 


2 


The  quantity  m  is  the  modulation  index  in 


that  it  is  the  ratio  of  the  rms  frequency  deviation  to  the  center 
frequency  of  the  modulation.  With  the  use  of  the  definition  of 

(A 

0 

Q  =  —  the  correlation  function  R  (t)  becomco 

<a).  V 


Before  attempting  to  evaluate  this  expression  it  will  be  helpful  to 
consider  limiting  cases  of  the  spectrum  S^(w)  as  m  becomes  very 
small  or  very  large. 


D.  LIMITING  CASES  FOR  FM  BV  BAND-PASS  NOISE 
1 .  Small- Modulation  Index 

Starting  with  the  expression  for  the  correlation 

function  of  v(t)  in  terms  of  S  ^  (w),  the  power  density  spectrum 

n 

of  the  modulating  noise,  we  can  obtain  an  expansion  of  the 
spectrum,  Sy(t*>),  that  is  appropriate  for  small  values  of  rrr. 
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by  the  S  (u)  under  consideration.  S  (u)  is  now  obtained  by  taking 
'^n  '' 

the  I'ouricr  transform  of  R^,{t) 


C«  2  2  r  r  -i. 

Sv(“)  =  “T  ®’‘P  1 


^2  2  r  Sv  ®v  1  1 

ni  ti  I  n  n _ _  4.  , . .  L 

L  (lo+w  (w-wl?~  J  J 


(W+Wg) 


(w-w  ) 


For  small  m  the  higher  ordered  terms  (tn^,  etc. )  can  be 

neglected:  thus  the  resulting  power  density  spectrum  of  v{t)  is 

seen  to  contain  a  carrier  spike  at  u  and  sidebands  about  the 

carrier  which  are  simply  ’•elated  to  the  modulating  noise  power 

spectrum  S  .  The  case  for  small  m  but  with  S  (w)  such  that 
n  '^n 

the  integral  1  does  not  converge  can  also  be  treated.  See 

Reference  5. 


2.  Large  Modulation  Index 


In  the  final  expression  for  Ry(T)  in  Section  IV-C 


let  a  =  and  w^t  -  t.  Then 


t  C.“  w_t 


exp  ^-m^j^l -exp(-a|t|)(cos  bt  +  ^  8inb|t|)|^. 
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It  is  clear  that  for  large  values  of  m  the  correlation  function 
becomes  very  small  as  1 1|  increases.  Therefore  a  good  ap¬ 
proximation  of  the  correlation  function  will  be  obtained  if  the 
terms  in  the  bracket  are  expanded  in  their  power  series  and  only 
the  lower  order  terms  are  retained.  Let  the  bracketed  term 
equal  f(t),  that  is, 

f(t)  r.  i  .,  exp(-a|  t|  )  (cos  bt  +  ^  sin  b|  t|  ) 

Expanding  the  exponential  and  the  cosine  and  sine  in  their  power 
series  gives 


,  3  .3, 

f(t)  =  1  -{l.a|t|  +  ^p -•••)(! 


=  1 


r ,  1  ,  2  ,2,  J.  , 

-  -  -jr  (a  -  b  )  t  + 


1  ,  2  .2.  J.  . 

=  7  (a  -  b  )  t  +  • '  • 


But  a^  =■  — and  b^  =  1  -  —^-<r  i  thus 


4Q‘ 


4Q 


f(t)  =  4 

2 

when  terms  jf  higher  order  than  t  are  neglected. 
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The  correlation  function  now  becomes 

-  2  _2  ? 

C  -  D  V  j  o  T 

«  <  .  0  I  n  1  2  2 , 

\{t)  =  -j-  cos  exp  ( - j-  ^  Uq  t  )  , 


Z  Z 

or,  if  (Aw)  =  D  where  Aw  is  the  mean  square  deviation  of 

the  carrier  w  , 
c 


r>  1  V  ^0  .  r  2  ' 

i.  J 


The  power  density  spectrum  is  now  obtained  by  use  of  the  Fourier 
cosine  transform 


S^(w)  =  -0-  J  exp  j-ilil^jcos  w^ 

-00 


T  COS  WT  dr  • 


Using  the  trigonometric  identity  for  the  product  of  cosines  gives 

.2_2 


C  2  z  Z 

Sy(“)  =  -^  J  e’tP  j  cos  (w  +  w^)t  dT 

-00 


*^0  r°°  r  (^w^)  1  ,  ,  J 

^  J  exp  -.1 - ^ —  fcos  (w  -  w^)t  dr 

-00 


2.  2 


s;(co)  +  s;(w)  . 
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If  Au  «  then  there  is  negligible  contribution  to  the  spectrum 
S  \<*))  by  S~  for  w  >  0,  and  similarly  S’(u)  is  very  nearly  equal 

J. 

to  Sy(w)  for  u  <  0.  It  suffices,  then,  to  consider  only  S^(«). 
From  Pierce's  Table  of  Integrals  No.  508  (6) 


s>) 


The  validity  of  this  result  is  easily  checked  by  noting  that 
Sy(u)  =  +  S'(w)  =  S'^(w)  +  S‘''(-w)  and 


Thus  the  power  in  the  carrier,  is  effectively  redistributed  into 
a  continuous  spectrum  by  the  process  of  frequency  modulation. 
Notice  that  the  spectrum  of  the  modulated  wave  v(t)  has  becor  e 
Gaussian  as  the  modulation  index  m  oo  and  thus,  has  the  same 
form  as  the  probability  density  function  of  the  modulating  voltage  v^{t). 
Middleton  (5)  has  shown  that,  under  quite  general  conditions,  the  power 
density  spectrum  of  a  frequency  modulated  wave  tends  to  assume  the 
shape  of  the  first  order  probability  density  function  of  the  modulating 

5 

voltage  when  the  modulation  index  goes  to  infinity  . 


See  also  References  7,  8,  and  9. 
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Now  that  the  cases  of  both  small  and  large  modulation 
indices  have  been  examined,  it  is  clear  that  the  power  density 
spectrum  of  FM  by  band-pass  noise  undergoes  a  rather  striking 
transition  as  the  modulation  index  is  increased  as  is  shown  in 
Figure  4.  The  question  to  be  considered  in  the  remainder  of  this 
chapter  is  the  spectral  shape  that  is  obtained  for  intermediate 
values  of  the  modulation  index,  m,  as  a  function  of  the  bandwidth 
of  the  modulating  noise,  or  equivalently,  the  filter  Q.  Ultimately 
it  is  desired  to  obtain  a  measure  of  the  deviation  from  the  limiting 
Gaussian  spectrum  that  will  result  for  various  values  of  m  and  Q. 

E.  DIGITAL  COMPUTATION  OF  FM  BY  BAND-PASS 
NOISE  POWER  DENSITY  SPECTRA 

In  order  to  provide  a  basis  for  comparison  with  power 
density  spectra  that  are  obtained  by  use  of  an  approximation 
technique,  the  Fourier  transform  of  R^(r)  was  performed  by 
digital  computation  for  selected  values  of  m  and  O.  Returning 
to  a  form  given  previously  on  page  48. 
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where 


f(t)  =  1  -  exp  {-a|  t|  )  (cos  bt  +  |  sin  b|  t|  )  , 

we  bav5  that 


c°°  CO 

(w)  =  J  R^(t)  exp  (-j(JT)  =  J  ^ 

-00  -CO  ®  ^ 


dt 

'0  '*'0 


r  expf-n'.  f'*-;]  exp  f  -j  ~ ]  dt  +  ^  f 

J  u  -I  L  <<^0  J  ^  J 


r  2  T  r  )t 

expj^-m  f(t)  j  exp[-j — j-?.  J  dt  , 

Again  if  the  assumption  that  »  Au  is  made  it  will  suffice  to 
consider  only  one  of  these  integrals.  The  second  integral,  which 
will  give  the  behavior  of  the  spectrum  around  w  =  u  and  li  = _ - 

C  '  Wjj  ' 

becomes 


Cq  .oo 


~  j  exp  |^-m^f(t)J  exp  (-j(5f)  dt  . 

“CO 

Note  that  this  integral  does  not  oxist  in  tho  stri;;t  oonso  becaudft 


lim  r  2,,  2 

t^O  ®’‘P[""’  =  cxp(-m‘‘)  ; 
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The  starting  point  in  this  analysis  is  again  the  correlation 
function  of  v(t)  which  is  reproduced  below  for  convenience: 


where 


with 


to  t 
c 


RJ^)  =  exp  j^-m^f(t)J  cos  ~ 


f(t)  =  1  -  exp  (-a  1 1  j  )  (cos  b 1 1 1  +  ^  sin  bhl)  , 


a  = 


■2tr  ' 


b  =  /I  - 


4Q‘ 


and  t  =  cjqT 


It  can  be  readily  seen  that  the  envelope  of  R  (-—  ),  that  is, 

C  2  -  2  n  “0 

—  exp  j^-m  f(t)J,  which-is  an  even  function,  has  successive 


maxima  at  the  points  1 1|  =  and  minima  at  ]  t|  = 
with  n  =  0, 1, 2,  •  •  • .  Attention  will  be  concentrated  on  the  behavior 
of  exp  j^-m^f(t)J  in  the  vicinity  of  the  maximum  points.  The  im¬ 
portance  of  these  points  in  determining  the  shape  of  the  resulting 
power  spectrum  becomes  obvious  when  one  observes  that  the 
envelope  function  decreases  lapidly  on  either  side  of  the  maximum 
points  even  for  relatively  modest  values  of  m,  for  example  on  the 
order  of  four. 
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To  determine  the  behavior  of  the  envelope  near  the  critical  points 
for  t  >  0  let  X  =  t  -  ;  then 


expjj-m^f(t)J  =  exp^-m^^  1  -  exp  (  )  exp  (-a|x|  ) 

l^cos  (b|  x|  +  2mr)  +  ^  sin  (b|x|  +  ?.mT)’J  I 


=  exp|^-m^  ji  -  exp(-^^^^)  exp(  -a|x|  )(cos  b|x|  +  -^  sinb|x|  )Jj. 

The  terms  exp(-a|xj),  cosbjxj,  and  sin  b|x|  are  expanded  in 
their  power  series  as  in  Section  IV-D-2  to  obtain 

2 

expj^-m^f(t)J  =  exp  |-m^  [^1  -  exp  (-  )  (1  -  +  '  '  •)]  j- 

-  exp|-m^[^l  -exp(-il^)jj.expj^-i^  exp(-.^^)x^J 

when  higher  order  terms  in  the  resulting  power  series  in  x  arc 
neglected.  Using  this  approximate  expression  for  the  envelope 
function  gives 


.  C«  w_t  CO 

=  T" 


^n?0  SW  . 
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where 


2  2_  2  .2 

r  ni  ,  2mraw.  Birnrii  f  .  ZTrnav/.  .  ^irn  »"] 

g(t)  =  expU-j-  exp( — g — Ht--g-)J+expj--2-  exp(-— g — ){t  +  )  | 


£  =  1  for  n  =  0  and  c  =2  for  n  =s  0 

n  n 


and 


For  cases  of;  interest  b  ^ - very  nearly  equal 

to  unity  since  values  of  Q  on  the  order  of  ten  or  more  will  be  con¬ 
sidered.  The  power  density  spectrum  is  found  to  be^ 


=  'n  -exp(-^)2j'X 


exp  ^^-— 2^  exp  (^)J  cos  Zimp 


where  p  =  - 2-  as  before.  Relatively  few  terms  of  this  power 

“o 

serir',  need  be  computed  to  determine  the  shape  of  the  power 
spectrum.  For  example,  the  curves  shown  in  Figure  6  were 
obtained  using  only  the  first  tliree  terms.  The  accuracy  o’  the 
approximation  is  readily  apparent  when  compared  with  the  digitally 
computed  curves  in  Figure  5. 


6 


See  Appendix  C. 
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As  a  smooth  spectrum  condition  is  approached  by 
increasing  m  with  Q  fixed,  the  major  contribution  to  the  devia¬ 
tion  of  the  spectrum  from  the  limiting  Gaussian  spectrum  is 
caused  by  the  first  maximum  in  the  envelope  of  the  correlation 
function  on  each  side  of  t  =  0,  that  is,  for  t  =  ±  ~  .  This 
corresponds  to  the  n  =  1  term  in  the  power  series  expression 

Also  the  maximum  deviation  for  the  Qaussian  limiting 
spectrum  occurs  at  (1=0.  Thus 

=  2exp{^)  exp{-m^[l  .exp(--J)]}  , 

where  AS^  ina.x  maximum  deviation  from  the  Gaussian 

spectrum. 

If  is  fixed  this  expression  provides  a  means 

of  determining  the  required  m  for  a  given  Q.  Two  experimentally 

determined  power  spectra  are  shov/n  in  Figure  7  illustrating  the 

appearance  of  the  spectrum  before  and  after  smoothing  takes  place. 

The  data  presented  in  figure  8  were  deter:nined  by  observing  the 

minimum  m  required  to  produce  an  arbitrarily  smooth  spectrum 

for  various  values  of  Q.  The  solid  line  was  computed  by  use  of 

the  relation  between  m  and  Q  for  AS^(3)  =  1  per  cent. 

v'^^'max  ^ 
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(a)  Power  Spectrum  Before  Smoothing 
m  =  4.0  fo  =  !2.l5mc  f,  =  0.6mc 
Horizontal  Scale  20mc/division 


(b)  Power  Spectrum  After  Smoothing 
m  =  6.5  fQ=l2.0mc  f|S0.6mc 
Horizontal  Scale  20mc/division 

FIGURE  7  FM  BY  BAND-PASS  NOISE  POWER  DENSITY 
SPECTRA  (EXPERIMENTAL) 
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FIGURE  8 


n  VERSUS  MODULATING  NOISI 

Q  FOR  A  SMOOTH  SPECTRUM  CONDITION 
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V.  FM  BY  BAND-PASS  NOISE  AND  A  SINUSOID 

In  this  chapter  the  approximation  technique  used  to  obtain 
the  power  density  spectrum  of  FM  by  band-pass  noise  will  be 
employed  to  study  the  case  of  FM  by  a  waveform  which  consists 
of  the  linear  sum  of  band-pass  Gaussian  noise  and  a  sinusoid. 
We  start  with  the  general  expression  for  the  FM  correlation 
function,  that  is 

^“c’’  ’ 


where  the  modulation  function  (j>(t)  is  now 

<t.(t)  =  D  J  "  ^m  ^ 


-00 

.t 


DE_ 


"  ^  I  cos  +  q) 


•00 


Thus,  the  modulation  function  is  seen  to  consist  of  a  term  due  to 
the  noise  identical  to  that  employed  in  the  FM  by  band-pass  noise 
only  case,  and  a  term  due  to  the  sinusoid.  Let 


and 


.t 

't'tt)  =  I  V  (t')  dt' 

i)  u 
-00 


r(t)  =  m  cos  (w^t  +  a)  , 

8  V 


DEm 

where  m  =  -  is  the  modulation  index  corresponding  to  the 

s  Wq 
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sinusoid  and  a  is  a  random  phase  angle  that  is  uniformly  distributed 
over  0  to  Uir.  Note  that  m^  is  defijied  in  the  usual  way  for  FM 
by  a  sinusoidt  that  is,  the  ratio  of  peak  deviation  and  modulating 
frequency.  The  correlation  function  can  now  be  written 

C  2 

C  * 

=  -2“  Rg  -[expUaj^T)  E  j^expOr^- jr^)J  E  [^exp  (j't'2-  j'l'i)]} 

with  the  last  step  being  Justified  by  the  fact  that  the  expected  value 
of  the  product  of  two  independent  random  variables  js  equal  to  the 
product  of  the  expected  values.  Note  that  the  shortened  notation 
r(t2)  =  r^i  etc.  has  been  employed. 

This  expression  for  the  correlation  function  is  the  same  as 
that  obtained  for  FM  by  band-pass  noise  only,  except  for  the  addi¬ 
tional  expected  value  term  due  to  the  sinusoid.  This  expected 
value  will  now  be  evaluated . 

E  fexp  =  Ejexp  (jr^)  exp  (-jr^jJ 

=  E  cos  (topt^f  a)  j  exp|^-jmg  cos  (“0*^1+  »)]}• 
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T’'e  exponentials  can  be  expanded  by  use  of  the  following  identity: 


00 


exp(jm  co8C)=  ^  e  J.  (m  )  cos  kO  , 


k=0 


k'“‘s' 


th 


where  is  the  k  Bessel  function  of  the  first  kind.  The  product 
is  written  as  a  double  sum  to  obtain 

Efexp  (jr2  -  jr^)J  = 

{  00  00  ,  *  ^  ^  ^ 

If  the  product  of  the  cosine  terms  is  expanded  to  obtain  terms  con¬ 
sisting  of  the  cosine  of  the  sum  of  the  arguments  and  the  difference 
of  the  arguments,  it  can  be  seen  that  the  bperation.of  carrying  out 
the  expected  value  integration  with  respect  to  the  random  phase  a 
will  eliminate  the  cosine  sum  term,  with  the  result  that 
Ej^exp  (jr2  -  reduces  to 

[(k-f)<^0t  +  ku>jjT  +(k-f)a]J. 


Further  it  is  seen  that  this  expression  yields  a  non-zero  value  only 
when  k  =  # . 
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Thus 


Using  the  results  obtained  for  FM  by  band-pass  noise  only,  we  find 
that  the  adding  of  a  sinusoidal  modulation  causes  the  resulting  power 
density  spectrum  to  become 


(n's)  cos[yp.k)], 


where  the  notation  of  Appendix  C  has  been  used. 

A  relation  between  the  modulation  index,  m,  which  corre¬ 
sponds  to  the  noise  modulation,  and  m^,  which  corresponds  to  the 
sinusoidal  modulation,  can  be  obtained  in  terms  of  the  ratio  of  the 


sinusoid  to  noise  average  power.  Let 


E  sin  (u>_t  +  a)  E  ^ 
m  '  0  '  r 


m 


V 

n  '  ' 


2  V 


D^E 

m 


2D^  V 


m 


2m 
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or 


m 


m  = 


The  power  density  spectrum  can  now  be  expressed  in  terms  of  the 
parameters,  m^,  p,  and  Q. 

2  .  „2 

E  J^(m  ) 
n  k'  s' 

,2 


•r^r4i~  ^ 


s  k=-CD  n=0 


p  { g )  expj-j:^  [l  -  exp  { -  j|> 


exp 


p(P-k)^ 


m 


exp{^) 


COS 


I^Zirn  (p-k)J 


where,  as  in  Chapter  IV,  the  assumption  has  been  made  that 
b  = 


1  '  «1. 


40“ 

This  expression,  although  rather  form’dable  at  first  glance, 

does  yield  an  insight  into  the  spectral  shape  that  is  obtained  from 

FM  by  band-pass  noise  and  a  sinusoid.  The  spectrum  is  seen  to  be 

made  up  of  the  sum  of  a  set  of  FM  by  band -pass  noise  spectra  which 

are  centered  at  integer  values  along  the  p  axis  with  each  having  a 

peak  value  corresponding  to  J,  ^(m  ).  The  degree  to  which  each 

X  s 

component  approaches  a  smooth  Gaussian  curve  will  depend  on  the 

m 

modulation  index,  m,  or  equivalently  -  .  Suppose  p  is  fixed; 

-i/^p 

then  the  value  of  m^  required  to  produce  smooth  Gaussian  shaped 
components  of  the  resulting  spectrum  can  be  determined  from  the 
results  of  Chapter  IV.  In  the  limit,  as  m^  becomes  large,  the 
spectrum  assumes  a  shape  given  by  the  probability  density  function 
of  the  sum  of  the  sinusoid  and  the  band-pass  noise. 
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VI.  SIMULTANEOUS  FM  BY  BAND-PASS  NOISE  AND  AM  BY 
LOW -PASS  NOISE 

The  waveform  to  be  considered  in  this  chapter  consists  of 


where 


and 


v(t)  =  C(t)  cos  l^w^t  +  <|>(t)J 


C{t)  =  Cq  +  v^(t) 


=  D  J  v^(t')  df 


-00 

The  voltage  i®  band-pass  noise  as  before,  while  v^{t)  consists 

of  noise  which  is  obtained  by  passing  white  noise  through  a  low-pass 
RC  filter. 

The  power  density  spectrum  of  the  low-pass  noise,  v  (t),  is 

a 


2  a>  V 

a  u>  +  w 


where  .  As  in  the  band-pass  noise  case  the  spectrum  has 

been  adjusted  so  that  a  constant  rms  voltage, ,  is  obtained 
independent  of  w^.  Note  that  co^  is  the  conventional  3  db  bandwidth 
of  liic  noise  processing  RC  filter.  The  correlation  function  of  v  (t) 

di 
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is  easily  obtained  by  taking;  the  inverse  Fourier  transform  of  (w) 

a 

to  get 


”  ''a'^  exp  (-10^1  t|  )  . 

a 

In  Chapter  III  the  correlation  function  of  simultaneous  FM  by 
band-pass  noise  and  AM  by  low-pass  noise  was  found  to  be 


Ry(T)  =  j  E|c(t^)  vJjL^)  cos  +  iKtj)  -  +(‘1)]} 


Substituting  for  the  amplitude  modulating  function  yields 


R„(,)  =  *  e([cJ  t  C„  v,(.,  )*  C„  v^C^l  +  v^(t2)  v^(.,)] 


+  t^t^)  - 


Using  the  fact  that  v  (t)  has  a  zero  mean  value,  and  the  assumption 

ct 

that  v^(l)  and  arc  independent  gives 


-Vt)  =  4{Co^  +  E[v^(t2)  v^{t^)]}  E  [cos  -  4>(ti)]}  . 

The  second  expected  value  is  again  identical  with  that  obtained  in 

Chapter  IV,  while  the  term  eTv  (t^)  simply  (t). 

U  J  a 
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Using  the  notation  o£  Section  IV-E,  we  get 


R  (1~\  -  ° 


jl  +11  exp  {-  -^--)1exp 


where  |j.  =  —  is  a  measure  of  the  percentage  amplitude  modulation. 

The  power  density  spectrum  of  simultaneous  FM  by  band-pass 

noise  and  AM  by  low-pass  noise  is  thus  seen  to  consist  of  two  terms: 

one  which  is  identical  with  that  of  FM  by  band-pass  noise  only;  and 

one  which  results  from  a  correlation  function  that  is  modified  by  the 

“a  1*^1 

factor  (A  exp  ( - ). 

a>Q 

The  effect  of  simultaneous  amplitude  modulation  can  be 

readily  estimated  by  considering  the  influence  of  the  factor 
U)  Jt] 

exp  ( — - - )  on  the  previous  derivation  of  the  spectrum  of  FM 

*^0 

by  band-pass  noise  only.  Recall  that  the  correlation  function  of 
FM  by  band-pass  noise  was  approximated  by  a  set  of  displaced 
Gaussian  functions  whose  amplitudes  were  adjusted  to  equal  the 
maxima  of  the  envelope  of  the  exact  correlation  function  at  the 

points  t  =  .  The  major  effect  of  the  factor  exp  ( — - ) 

is  to  reduce  the  peak  amplitudes  of  the  approximating  functions. 

Thus  the  corresponding  effect  on  the  spectrum  becomes  one  of 
reducing  the  magnitude  of  the  degree  to  which  the  r.ower  density 
spectrum  deviates  from  the  limiting  Gaussian  shape.  With  these 
considerations,  the  results  of  Appendix  C  for  FM  by  baod-pass 
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noise  can  be  readily  extended  to  the  case  of  simultaneous  FM  by 
band-pass  noise  and  AM  by  low-pass  noise  to  obtain 
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VII,  EXPERIMENTAL  INVESTIGATION  OF  POWER  DENSITY 

SPECTRA 

The  system  used  to  generate  the  vai'ious  random  modulated 
waveforms  discussed  in  thin  report  is  shown  in  Figure  9>  Spectrum 
analysis  was  accomplished  using  a  modified  AN/APR-9  receiver, 
which  is  essentially  a  double  IF  frequency  superheterodyne  receiver 
with  RF  preselection  employed  for  image  rejection.  The  local 
oscillator  klystron  frequency  is  linearly  swept  between  adjustable 
frequency  limits  by  a  mechanical  drive  unit.,  A  display  of  the  power 
density  spectrum  of  the  receiver  input  waveform,  such  as  that 
shown  in  Figure  7,  was  obtained  by  use  of  an  auxilliary  oscilloscope 
with  a  long  persistence  cathode  ray  tube.  A  voltage  proportional  to 
the  receiver  local  oscillator  frequency  was  applied  to  the  oscillo¬ 
scope  x-axis  input:  while  the  output  of  the  receiver's  narrow  band 
second  IF,  after  suitable  detection  and  integration,  was  app'ied  to 
the  y-axis  inpiJt  of  the  oscilloscope.  This  method  of  rapid  power 
spectrum  moasuremeni.  rriakes  feasible  an  investigation  of  the  effect 
of  the  various  modulation  parameters  that  would  be  difficult,  if  not 
impossible,  if  the  conventional  point-to-point  power  density  spectrum 
measurement  technique  were  used. 

An  accvirate  determination  of  the  frequency  sensitivity  of  the 
FM  source  (BWO)  in  radians/volt  (denoted  by  D  in  tins  report)  was 
obtained  by  applying  a  sinusoidal  modulating  volt.age  of  sufficient 
amplitude  to  cause  the  center  frequency  spectral  line  to  he  r.evo. 


FIGURE  9  BLOCK  DIAGRAM  OF  EXPERIMENTAL  EQUIPMENT 
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2 

Since  the  center  frequency  spectral  line  is  given  by  Jq  {^g)i  the 
first  null  of  this  spectral  line  will  occur  at  a  value  of  corre¬ 
sponding  to  the  first  zero  of  1^(111^).  The  peak  magnitude,  of 

the  modulating  sinusoid  that  is  required  to  produce  the  first  null  of 
the  center  frequency  spectral  line  is  measured,  and  D  is  then  com- 
puted  from  the  relation  D  =  — ^  ,  where  f^  =  -j-  is  the  frequency 

of  the  modulation.  A  similar  technique  using  other  spectral  lines 

2 

having  magnitudes  of  Jj^  used  to  determine  the  value 

of  D  corresponding  to  various  modulation  voltage  magnitudes  and 

thereby  determine  the  modulation  linearity  of  the  FM  source.  Once 

— 2 — 

D  has  been  determined,  the  mean  square  noise  voltage  v^  (t),  re¬ 
quired  to  produce  a  desired  noise  modulation  index,  m,  is  computed 
from  the  relation 


2  2 
u>0  m 


''n  -  ' 


and  v^(t)  in  adjusted  to  obtain  the  required  w.ith  the  aid  of  a 


true  rms  meter. 
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APPENDIX  A 


If  a  band-pass  circuit,  whose  transfer  function  is 


H(8)  =  --2 - - - 2- 

S  +  <J^S  +  Wq 

is  used  to  filter  white  noise,  that  is,  noise  having  a  constant  power 
density  spectrum,  then  the  resulting  output  voltage,  has  a 

power  density  spectrum  given  by 


S„  («)  = 


Hq^  N 


.1  ,  2  2,  .  2  2 

(Wq  -  W  )  +  <»> 

where  N  is  the  constant  power  density  level  of  the  input  noise.  The 
mean  square  output  voltage  is 


Hq'N  oo 


2  j 

u  du> 


"ZF 


•'  ,2  Z,  ^  2  2 

-00  (Wq  -  W  )  +  u> 


This  integral  can  be  readily  evaluated  by  contour  integration. 
By  analytic  continuation  the  integral  becomes 


I 

C 


.2,2,^  22 
(W-  +  7.  )  +0)^  z 


dz 


where  C  is  a  closed  contour  along  the  real  line  -oo  <  w  <  oo  and 

encircling  the  tipper  half  plane.  The  integrand  has  two  poles  that 

/  2 

lie  within  the  contour,  one  at  /,  =  “''./“o  ”  ~T~  ’*'  •*  "Z~ 
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another  at  z  =  ,  The  integral  can  be  written 


z^  dz 


C  )(z+Z|)(zi+z^*) 


=  2irj  l^f 


Sum  of  Residues  at  and  -z^  J 


2  .  ^ 

-  ’.iF  ,  '■“*  ’  1 

-  - 5f - 5“  + - r - s - ar - s - 

)  {-*1  -z^^)(-z^*+z^)  J 


=  If  j 


(Zj+iSl  ) 


*~1  '  ~1  '  „  _  Tfj  -  _  irj  ^ 

2  >!<2  ‘  z]  -  z/  "  "  “l 

1  "  ^  ^ 


Note:  The  symbol  Im  denotes  "imaginary  part  of,  " 

The  integral  over  the  upper  half  z-plano  circle  can  be  easily 
shown  to  be  zero  in  view  of  the  fact  that  the  degree  of  the  denomina¬ 
tor  of  the  integrand  is  greater,  by  two,  than  the  degree  of  the 
numerator.  Thus  the  result 


-00 


2  , 

(0  du) 


,2  Z,  ^  2  2 

(Wp  -  10  )  +  (0 


is  obtained. 
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APPENDIX  B 

In  order  to  obtain  the  correlation  function  of  FM  by  band-pass 
noise  an  integral  of  the  form 

I(w)  =  r”  _ 

*'  ,2  2,  .  2  2 

-00  (Wq  -  O)  )  +  U)^  U) 

must  be  evaluated.  The  integral  I{w)  can  be  written  as  the  sum  of 
two  integrals,  IA<^)  and  l,(t>)),  with 


and 


i(“i  =  I 


du 

T 


,2  2,  2  2 

-00  (Wp  -  U)  )  +  w 


I 


,(»)  =  J 


-cos  wT  du 


2"'  ■  J  ,  2  2,2  TT 

-00  (Wp  -  w  )  +  w 


These  integrals  are  evaluated  by  means  of  contour  integration. 
By  analytic  continuation,  the  first  integral,  I^(w),  becomes 


i'“>  -  J 


dz 

T“ 


^  ,  2  2,“  2  2 

C  (Wp  -  z  )  +  z 

where  the  contour  is  taken  along  the  real  line  -oo  <  to  <  oo  and  the 
upper  half  plane  with  the  integral  around  the  upper  half  plane  being 
zero  as  was  the  case  in  Appendix  A.  Note  that  the  poles  of  the 
integrand  are  also  the  same  as  those  of  the  integral  in  Appendix  A. 
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I^(u)  =  2irj  I^Sum  of  Residues  at  and  -z^  *'J 


,  2 
<*r*i 


I  1  X  1  » 

)  ^  ""l  “i“l 


z,z’!‘{z.-z*) 


"l“0 


The  integral  I,{w}  is  evaluated  as  follows,  where  the  symbols 
Re  and  Im  denote  "real  part  of"  and  "imaginary  part  of" 
respectively: 


-I2H  = 


r- 

-00  (tij 


COB  (a)T  dto 

2  z7~~Tz 

O'"  )  +"i" 


exp  (jzlr!  )  dz 

Tz  zJTzz 

(Wq-z  )  +  z 


=  Re|- 


JLL 


fexpljz^lrl)  exp  (-j  z*1t|)  1 


2 

*1  “  *1 


} 


=  Re 


( - ^ - 2-  r  **exp(jzjT|)  +  7,^  exp(-jz*|T!)1  i 

*.  I_  I  Z#_2  ^  J 


Zj  (Zj-Z^  ) 


=  Re 


— -7-31 - exp  R-I  rRe(z  )coiRe(zJT!) 

l2|zJ^Re(z^)IJz^)  L  rn  1  JL  1  1 
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APPENDIX  C 


In  Chapter  IV  the  correlation  function  of  FM  by  band-pass 
noise  is  approximated  using  the  following  expression; 


W  t  GO 


J^r.^exp[-,n^[l-exp(-i’3^)1}.  g(t)  , 


where 


8(.|  =  e.p[-^exp(- + 


Let 


and 


then 


Also  let 


Z  2  ,  Zxna . 

*"  p  nr“^ 


.  Zirn 
‘n  -  "T"  ' 


g(t) 


Z  ^  Z 

=  exp  j^-  (t  ••  t^)^j  +  j^-  (t  + 


Kn  =  exp  /-m^  |^1  -  exp  (- 


Now  can  be  written 


w  t  00 


R  ^  COB  — ^  y*  e  K  g(t) 

v  “0  '  ^  <-0  n=0  "  ^ 


The  power  density  spectrum,  Sy(w),  can  now  be  obtained  by  taking 
the  Fourier  transform  of  R  (7^)-  Thus 


S  (w)  =  f  R  {  — )  exp  {■— )  — 

v'  '  J  v'wq'  ‘^0  0 

-00 

“  *  t 

Upon  writing  the  cos  -2-  term  contained  in  R  ( — )  in  exponential 
'•*0  0 

form,  Sy(w)  becomes 

00  (•  00  (w-w  )  -,  «<»  r  n 

J. ‘'„{J  g(t)e«p[-l  — t/l+J  gWexpCj— t] 

“  ^rv\  -•nr^ 


As  before  (Section  IV -E)  we  assume  that  »  oj^,  which  leads  to 
the  conclusion  that  the  power  density  spectrum  appears  centered 
about  the  frequencies  and  -w^.  Since  8^(0^)  is  an  even  function, 
attention  need  only  be  focused  on  the  spectral  behavior  around 
w  =  which  is  obtained  by  using  only  the  first  integral  in  the 
expression  for  Sy(w).  Thus  S^(P),  which  is  S^(w)  expanded  about 


Oi)  =  u  ,  with  0 
c 


C  ^  CO 

0  r> 


“>'-4-  X'nKn  J  8W«xp(-iPt|<l, 

0  n=0  " 


An  asymptotic  power  series  approximation  for  the  power  density 
spectrum  of  FM  by  band-pass  noise  is  thus  obtained  upon  evaluation 
of  the  integral 


g(t)  exp  (-jpt)  dt 
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Substituting  for  g(t)  leads  to  the  two  integrals 
2  CO 

/”  [‘  "T  J  [■  "T  jP*]  ' 


If  the  first  integral  is  evaluated,  the  second  integral  can  be  readily 
obtained  by  simply  changing  the  sign  of  t^.  The  first  integral  is 
seen  to  be 


J^exp|^--^(t^+ 2(t^+  t+ t^^J  dt  , 


which  upon  completion  of  the  square  becomes 


®’‘P  [■  ■~7  ■*■  j  *n^  J  ^’‘P  ■  -r-  ^  • 


u  =  t  +  t^  + 


then  this  integral  becomes 


r'®  r  %  2]  ^  -/I 

J  exp|^--2-  u  j  du  . 


This  last  integral  uses  the  well-known  result  of  the  integral  of  the 

2  1 

Gaussian  first  order  probability  density  function,  where  o-  =  — ^ 

a 

is  the  variance.  ” 
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rhiis  the  integral 


^oo 


I  g(t)  exp  (-jpt)  dt  = 


-00 


•J' Ztt 


^  )  [e^P  {j  P>  +  exp  ( -j  p)J  = 


2  yi7  ,  0^  . 

- -  exp{-.-K^5  cos  t  8 

n  2o  ^  ” 

n 

and  the  power  series  approxinnation  for  S^(P)  is 

®>|  =  iSj  •/!  J„  ^  »«p  ( -  t„  p 


"0 

2 


Jq  ""n  «’‘P(^)  exp|-m^|l-exp(.^)JJ 
exp(Ig2:)j  cos^  p 

Notice  that  the  zeroth  order  term  of  this  expression  is 

/r  ,  P^  . 

TSUV^  exp(-.P  )  , 

Zm 

which  is  the  same  as  that  obtained  in  Section  IV-D-2  for  the  limiting 
case  as  the  modulation  index  m  becomes  very  large.  Thus  the 
higher  ordered  terms  in  the  expression  for  S^(p)  represent  devia¬ 
tions  of  the  spectrum  from  the  Gaussian  limit  for  iafermediate 
values  of  m. 
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